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We study the kinematics of dark matter using the massive cluster of galaxies 1E0657-56. The 
velocity of the "bullet" subcluster has been measured by X-ray emission from the shock front, and 
the masses and separation of the main and sub-clusters have been measured by gravitational lensing. 
The velocity with gravity alone is calculated in a variety of models of the initial conditions, mass 
distribution and accretion history; it is much higher than expected, by at least 2.4a. The probability 
of so large a subcluster velocity in cosmological simulations estimated from the Hayashi- White fit 
is < 10" 7 . A long range force with strength ~ 0.4 - 1.2 times that of gravity would provide the 
needed additional acceleration. 



Determining whether or not dark matter (DM) has 
long range non-gravitational interactions will help decide 
whether it is just a weakly interacting particle which has 
not yet been observed in the laboratory or a manifes- 
tation of an entirely new sector of nature. Generically, 
extensions of the standard model such as supersymme- 
try and string theory imply the existence of scalar fields 
("moduli") which are massless at all orders of pertur- 
bation theory and naturally lead to long-range inter- 
actions called "fifth forces" for historical reasons. To 
extraordinary accuracy, gravity is the only long-range 
force affecting baryonic matter but constraints on 
non-gravitational forces within the dark sector are much 
weakcr[2] . A fifth force may be of a Yukawa form, approx- 
imately l/r 2 near the source and exponentially damped 
beyond some range, rs, which is most trivially a constant 
or may grow with the expansion of the Universe Q. In 
chameleon models [H the new interaction is damped by 
sufficiently strong concentrations of matter, so it is im- 
portant to test for the existence of a fifth force on many 
different scales, since it may either appear or disappear 
at large scales. In this work we examine the constraints 
on dark matter self-interactions on the Mpc scale, which 
can be obtained using a very special system, the colliding 
galaxy clusters 1E0657-56. Thanks to the proximity and 
fortuitous geometry of this system, it is possible to mea- 
sure the masses of the main cluster and the subcluster, 
and their relative velocity, and thus calculate whether the 
acceleration implied is consistent with the force between 
them being purely gravitational. 

An attractive non-gravitational force between DM con- 
centrations is not only well-motivated theoretically, it 
may resolve some discomforts with conventional ACDM. 
If the new force is long range compared to the scale of 
structure formation, it effectively modifies the strength 
of the gravitational interaction for the dark matter to be 
(l+P)MiM2 Gn and thus would accelerate structure for- 
mation. • The number of superclusters observed in SDSS 
data appears to be an order of magnitude larger than 
predicted by ACDM simulations @; accelerated structure 
formation would reduce this discrepancy. • As noted in 
Q , a fifth force would tend to clear out the voids; ref. @ 
confirms this in a simulation. This may improve agree- 



ment with ACDM [7], although the existence of a dis- 
crepancy is not certain^. • A variety of observations, 
for instance the lack of evidence in the Milky Way for 
a major merger, is hard to reconcile with the amount 
of accretion predicted in ACDM. Accelerated structure 
formation reduces late-time accretion, simply because it 
leaves less to be accreted later @- • The number of satel- 
lites in a galaxy such as the Milky Way is predicted to 
be an order of magnitude larger than is observed. This 
"substructure problem" is ameliorated by a 5th force, by 
reducing the stellar content of dwarf galaxies and mak- 
ing them harder to find. This is because baryons - not 
feeling the 5th force - are relatively less-bound to dark 
matter concentrations than in conventional theory, re- 
ducing the amount of bound gas and lowering the star 
formation rate in dwarf galaxies, and increasing the tidal 
loss of the stars that do form. 

Structure in the Universe evolves by merger and accre- 
tion, so many large clusters of galaxies are observed to be 
disrupted due to a recent major merger. A combination 
of factors makes 1E0657-56 unique: 

• The bullet subcluster has "just" fallen through the 
main cluster for the first time, so its baryonic components 
- hot gas and galaxies - are still substantially intact. 

• The subcluster and main cluster are both exceptionally 
massive, their relative velocity is very high, the system is 
nearby from a cosmological perspective (z — 0.296), and 
the trajectory of the subcluster is nearly exactly in the 
plane of the sky. Therefore, the supersonic shock front 
between the subcluster and main cluster gas is clearly 
visible in the X-ray images and the geometry is simple 
and clean enough to give an excellent determination of 
the relative velocity [lCj| . 

Clowe, Markevitch and collaborators have used 
1E0657-56 to exclude MOND-type alternatives to 
DMpj], [12], place limits on short-range interactions of 
DM particles a nd test models of X-ray production in 
massive clusters TO]. (See these references for a survey of 
earlier literature on 1E0657-56.) With colleagues, they 
have also undertaken a massive campaign to improve the 
X-ray observations and gravitational lensing constraints. 
A new 500 ks Chandra observation of the discontinuity 
in the density of gas across the prominent bow shock 
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preceding the gas bullet allows precise determination of 
the shock Mach number and according to [l0( implies 
a shock (and bullet) velocity 4740i"55 l kms~ 1 . The er- 
rors quoted in jToj were symmetrized and the central 
velocity rounded down(M. Markevitch, private commu- 
nication); using more accurate values results in a veloc- 
ity of 4740+55° kms -1 which we adopt here. Combining 
weak and strong lensing gives a much improved mass 
distribution in the inner 500 kpc region[l4| and leads to 
a separation between main cluster and subcluster mass 
peaks of 720 ± 25 kpc. A much larger weak lensing 
field (34' x 34') allows the mass distribution to be fol- 
lowed to larger distances. Applying the same weak lens- 
ing analysis described in the new large-field weak 
lensing data leads to (D. Clowe, private communication) 
^200 = 2136 kpc, c = 1.94 and r2oo = 995 kpc, c = 
7.12, for the main and subclusters respectively, taking 
them to have spherically symmetric NFW profiles with 
h = 0.7, Qm = 0.3. These parameters imply masses of 
A/ 200 = 1.5 x 10 15 M SU „ and Af 200 = 1.5 x 10 14 M stm for 
the main and subclusters respectively, and virial veloc- 
ity V200 = 1740 km s^ 1 for the main cluster. A King 
profile fits the weak lensing data slightly better, with pa- 
rameters po = 2.16 x 10 6 M sun kpc~ 3 , r c — 264 kpc and 
po = 1-01 x 10 7 A/ su „kpc~ 3 , r c = 78 kpc. More general 
mass distributions for the main cluster were also tried; 
they lead to lower total mass estimates. See 15| for a dis- 
cussion of systematic uncertainties in weak lensing fits. 

Given the data above, we determine the most proba- 
ble "fiducial" model and use it to calculate the expected 
velocity of the bullet subcluster. Variants on the fidu- 
cial model give a sense of the uncertainty in the esti- 
mate. If the main cluster's mass distribution were static 
and azimuthally symmetric about the direction of motion 
of the subcluster, the gravitational accelerations from - 
720 to +720 kpc would cancel. This trip takes less than 
about 0.4 Gyr, during which time mass accretion is a 
small effect, so our calculation is relatively insensitive to 
the mass distribution in the central region. This is for- 
tunate, because the central mass distribution is likely to 
be complicated due to this and earlier mergers, and may 
not be well described by an NFW profile. In the region 
between the central core and the virial radius, a spher- 
ical NFW density distribution fits simulations and data 
well on average, so we take our fiducial mass distribution 
for the main cluster to be a spherical NFW mass profile 
embedded in an otherwise homogeneous Universe of den- 
sity /?m,o(1 + z ) 3 > with parameters fixed by the Clowe fit 
quoted above. 

The mass distribution of the main cluster evolves in 
time due to mergers and accretion so we need the mass 
distribution not only at z = 0.296, but at earlier times 
as well. The mass accretion history (MAH) of galaxy 
clusters has been extensively studied [lg]. In most cases it 
is well-represented by the function presented in Wechsler 
et al [l6| and we adopt this for the fiducial case. To gauge 
the degree of dispersion we followed the actual MAHs 
of the 12 clusters with mass > 10 15 Af su „, in a recent 



simulation of ~ 80,000 galaxy clusters [17]. In the most 
extreme case, the velocity using the actual MAH was 10% 
higher than obtained using the mean MAH. 

The cumulative effect of the multiple small gravita- 
tional deflections a test body experiences when pass- 
ing through an ensemble of point masses with a dis- 
tribution of velocities is known as dynamical friction. 
At each position along the trajectory of the subcluster, 
we approximate the deceleration due to dynamical fric- 
tion by the Chandrasekar formula (see, e.g., Binney and 
Tremaine[l8j) integrated over impact parameters in the 
range relict to 6 max and assuming a locally Maxwellian 
velocity distribution. Here rbuiiet is a characteristic size 
of the bullet which we take to be r2oo/c = 140 kpc and 

6 max = (yj^Qt^j is the scale over which the local density 
changes by of order one; we checked that results are only 
weakly sensitive to these choices. 

Note that the initial time and position are not inde- 
pendent because the sub-cluster is constrained to be on 
its first exit through the center and to have reached a 
radius 720 kpc at z = 0.296. Hence the predicted final 
velocity is insensitive to the distance at which infall is 
presumed to start, because the acceleration is negligible 
until the excess mass in the main cluster above the cos- 
mic mean background density becomes significant. The 
rms peculiar velocity of galaxy clusters is 293 ± 28 km 
s _1 , with < 5 % probability of a velocity greater than 
600 km s _1 [l9[ . Therefore we conservatively take the 
fiducial value of the initial infall velocity to be 300 km 
s _1 , and also consider and 600 km s _1 . 

The fiducial model predicts a much smaller velocity 
than the 4740 kms -1 observed: 2950tgo° kms -1 , where 
the uncertainty range reflects the 0-600 kms 1 range of 
initial velocities. To produce the observed final veloc- 
ity, an absurd initial velocity of 3135 km s _1 would be 
needed. We varied the parameters and assumptions of 
the model to see how much the predicted velocity can be 
increased; Table U summarizes our results. We consid- 
ered a King profile and an NFW profile truncated at the 
virial radius; these give still lower subcluster velocities, 
the larger being 2330 kms -1 for the King profile. An 
N-body simulation would do a better job on dynamical 
friction, but as shown by the cases with no dynamical 
friction, improving the treatment should not substan- 
tially affect the conclusions since dynamical friction is 
responsible for only a 325 kms -1 decrease in the pre- 
dicted velocity. To illustrate the sensitivity to NFW pa- 
rameters, we fit the surface mass density within 500 kpc 
from the combined weak and strong lensing analysis of 
Bradac et al[l4j]. This is a promising approach but less 
secure for our application since it is sensitive to prob- 
able deviations from the NFW profile in the disrupted 
inner region, which is inessential for us. Furthermore, 
it requires an extrapolation relying on a strictly NFW 
profile, to the larger distances relevant to our analysis. 
With these parameters we obtain a larger velocity, but 
still much lower than observed. The subcluster velocity 
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is largest for a prolate NFW profile whose major axis is 
aligned with the direction of the subcluster. To maxi- 
mize the effect, we used c/a = 0.65, corresponding to 
the most elongated example among the very massive ha- 
los of [2(1; we re-fit the Clowe weak lensing profile to 
obtain 1E0657-56 parameters under this analysis. Even 
this most extreme case gives a velocity 3200 km s" 1 , much 
below that observed. Reducing the separation between 
bullet and center of main cluster by la to 695 kpc has 
a negligible effect, only increasing the predicted velocity 
by 35 km s _1 . 



TABLE I: Expected velocity of the bullet subcluster under 
various assumptions; cluster mass density model; notes. 



v S ub,{ (km/s) 


Profile 


Notes 


2950 


NFW, std 


c = 1.94, r-200 = 2.136 


2840 


NFW, std 


v in = 


3080 


NFW, std 


v in = 600 


2985 


NFW, std 


rbuiiet = 695 kpc 


3275 


NFW, std 


no dyn fric 


2330 


King 


no MAH 


2785 


King 


no MAH, no dyn fric 


3200 


NFW, prolate 


c = 2.16, r 2 oo = 2.235 


3530 


NFW, Bradac 


c = 5.22, raoo = 2.235 


3435 


prolate + unseen mass 


conspiracy model 



Could matter not included in the NFW profile, at dis- 
tances beyond the 3 Mpc of the weak lensing observa- 
tion, possibly exert a sufficiently strong additional force 
to account for the observed velocity of the subcluster? In 
the ROSAT image, 1E0657-56 seems isolated and there 
is no indication of other galaxy clusters in an extended 
filamentary structure. Given its sensitivity limit, the 
ROSAT observation excludes a cluster more massive than 
about one-tenth that of 1E0657-56 within about 10 Mpc 
(A. Vikhlinin, private communication). Other variations 
in the assumptions only increase the discrepancy with 
the prediction of gravity alone. If the motion of the clus- 
ter is not precisely in the plane of the sky, or if it does 
not pass directly through the center of the main cluster, 
the subcluster velocity extracted from the shock front 
discontinuity increases because the actual distance from 
the center of the cluster to the leading edge of the shock 
is then higher so the ambient temperature and density 
of the pre-shock gas is lower and the density contrast is 
greater. (This is the reason for the asymmetric velocity 
errors; we thank M. Markevitch and A. Vikhlinin for dis- 
cussions of this point.) Moreover, the subcluster would 
be further away from the center of the main cluster, so 
would have been deccelerated more. If the trajectory 
does not pass through the center of the main cluster, a 
given observed velocity implies a longer travel time and 
at any given radial position along the trajectory the age 
and mass of the main cluster are lower. 



With the above in mind, we create a "conspiracy" 
model which combines those features that enhance the 
predicted velocity, at their 95% confidence limits - ex- 
treme prolate mass distribution aligned in the direction 
of motion, infall velocity of 600 km s _1 , and a static 
mass 1.5 x 10 14 M sun located 3 Mpc away from the center 
of the main cluster in the direction of motion of the sub- 
cluster - and we reduce the observed separation by la 
to 695 kpc. This conspiracy model predicts an observed 
velocity of 3435 km s _1 . 

If the anomalously large observed velocity of the bul- 
let subcluster is the result of a new Yukawa-like interac- 
tion in the dark sector whose range is larger than several 
Mpc, its effect on the bullet dynamics is to multiply the 
gravitational acceleration by the factor (1 + (3). For the 
fiducial case, the value of j3 which would bring the pre- 
dicted velocity into agreement with the central value of 
the observed velocity for rg = oo is (3 = 1.2, and for the 
conspiracy case is (3 = 0.8; reducing the final velocity by 
la to 4190 km s _1 , (3 = 0.4 is needed in the conspiracy 
model. These values are generally consistent with the 
limits on a fifth force found in [2j which were only rough 
estimates due to the much more primitive state of obser- 
vational astrophysics at that time. Just as the dynamics 
of 1E0657-56 is insensitive to the mass distribution inside 
~720 kpc, it is likewise insensitive to a fifth force in this 
region and would be incapable of discriminating between 
a simple Yukawa interaction and a chameleon interaction 
which is suppressed in the center of the cluster. 

Hayashi and White in 21] (HW) have recently taken 
a complementary approach to the direct dynamical anal- 
ysis presented above, investigating how frequently sub- 
clusters in the Millenium Run simulation[5] have com- 
parable characteristics to the bullet. Using the bullet 
data prior to 2005, HW find the likelihood of a bullet- 
like subcluster to be about 1 in 500. The new obser- 
vations reduce this because fbuiiet/V^oo increases from 
1.9 adopted by HW to 4740/1740 = 2.7. HW fit the 
distribution of velocities of most massive sub-clusters di- 
vided by the virial velocity of their host, T4 u b/V2oo, to 
obtain: log ^gf) = _ (Kh^oo)^ with Vw% = L55 

and a — 3.3 at z = 0.28. Thus increasing the ve- 
locity ratio from 1.9 to 2.7 reduces the likelihood by a 
factor 4 x 10~ 5 so the HW likelihood estimate becomes 
0.8 x 10~ 7 . This remains small even when measurement 
errors are included: reducing ^bullet by la to 4190 km 
s _1 and increasing M 2 oo by 20% (a very conservative es- 
timate of the uncertainty, obtained by augmenting the 
quoted 16% error on the first small field, low statistics 
weak lensing mass determination [11] by 25% for system- 
atic errors) gives Vbuiict/V2oo = 2.2 for a 1.4 x 10~ 4 like- 
lihood. These are only estimates because the observed 
ratio fbuiiet /V200 is higher than any observed in the sim- 
ulation so the likelihood values rely on the HW extrap- 
olation formula. Nonetheless, they make clear that the 
bullet is rare enough it should arouse suspicion and de- 
serves further study. 
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Other manifestations of a non-gravitational force be- 
tween dark matter concentrations are often masked by 
our lack of independent information on the underlying 
DM system. An important and clean test would be to 
stack SDSS clusters to obtain their total mass distribu- 
tion from weak lensing, and to map the peculiar velocity 
distribution for the same stack of clusters. If there is a 
fifth force, the maximum excursion in the velocity dis- 
tribution will be larger than expected (22j, although a 
simulation would be needed to interpret the results given 
complications such as dynamical friction, stripping, etc. 
On a sub-galactic scale, Kesden and Kamionkowski have 
pointed out that a difference in acceleration between DM 
and stars would change the distribution of stars in the 
tidal tails of the Sagittarius dwarf galaxy [23| . They argue 
that the method should be sensitive to /3 > 0.1, but until 
simulations with gravity alone can consistently describe 
all observations it is premature to claim any limits on (3. 
In particular, it is intriguing that the difficulty of recon- 
ciling the observed line-of-sight velocities of the stars in 
the leading stream, which implies a prolate Dark Mat- 
ter halo, with the precession of the debris orbital plane, 
which implies an oblate haloJUj], might be resolved by an 
additional attractive force. 

We have shown that it is very difficult to reconcile the 
reported velocity of the bullet subcluster with the ob- 
served matter distribution which accelerated it. With 
a downward error of 550 km s -1 , the subcluster ve- 
locity is nominally 3.25a higher than the fiducial pre- 
diction of 2950 km s^ 1 . A conspiracy model, requiring 
95% CL extremes of initial infall velocity, elongation ratio 
of a prolate mass distribution, a hypothetical additional 
mass concentration at the estimated maximum level al- 
lowed by weak lensing and ROSAT, and a la smaller 
separation, predicts a velocity of 3435 km s _1 , still 
falling short of observation by 2Aa. Clearly, further ef- 
fort to reduce uncertainties in all relevant quantities is 
warranted. Reducing reliance on simple parameteriza- 
tions of the mass distribution would be particularly help- 
ful. Adding more red-shifts for the arcs and extending 
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observed dynamics of the bullet subcluster is substan- 
tiated by refined observations and analysis, and con- 
firmed in other systems, it would imply the existence 
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sector. This would have profound implications for parti- 
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range of the force, how it varies with the scale factor of 
the Universe, and whether the Yukawa description is ap- 
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simple exchange of an ultralight scalar particle, which are 
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that dark matter will be difficult or impossible to detect 
in laboratory experiments, because otherwise quantum 
corrections would give rise to an unobserved equivalence- 
principle- violating interaction of ordinary matter; details 
will be reported elsewhere. 
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